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We study the magnetic resonant x-ray scattering (RXS) spectra around the K edge of Cu in 
KCUF3 on the basis of an ab initio calculation. We use the full-potential linearlized augmented plane 
wave method in the LDA+C7 scheme, and introduce the lattice distortion as inputs of the calculation. 
We obtain finite intensity on magnetic superlattice spots, about three orders of magnitude smaller 
than on orbital superlattice spots, by taking account of the spin-orbit interaction (SOI). No intensity 
appears without the SOI, indicating that the intensity arises not from the spin polarization but 
from the orbital polarization in 4p states. The present calculation reproduces well the experimental 
spectra as functions of photon energy and of azimuthal angle. We also calculate the RXS intensity 
on orbital superlattice spots. It is found that the intensity increases with increasing Jahn- Teller 
distortion. The spectra remain nearly the same in the nonmagnetic state given by the simple LDA, 
in which the orbital polarization in the 3d states is much smaller. This strongly suggests that the 
intensity on orbital spots is mainly controlled by the lattice distortion, not by the 3d orbital order 
itself. 



I. INTRODUCTION 

The resonant x-ray scattering (RXS) has attracted 
much attention, since the orbital order can be directly 
probed by the RXS signals on the orbital superlattice 
spots. Several RXS experiments ^jse-befiii carried out 
on typical perovskite compounds ,EmfjuratJ0 in which an 
orbitally ordered state may be stabilized below a critical 
temperatureE'Ej For the K edge in transition-metal com- 
pounds, RXS intensities arise from the modulation of 4p 
states in the intermediate states of the resonant process. 
Since Ap states are not the states of orbital ordering, this 
causes complications on the interpretation of RXS inten- 
sities. 

A mechanism was proposed that the modulation of 
Ap states comes from orbitally polarized 3d states on 
the Mn 3+ ion through the anisotronkjpterms of the 3d- 
4p Coulomb interaction in LaMn03,E3liil and thereby it 
was argued that the RXS intensities on the superlattice 
spots are a direct reflection of orbital order. On the 
other hand,, sub sequent studies based on band structure 
calculational3ll3ll3 have revealed that 4p states of Mn 3+ 
ion so extend in space that they are considerably modified 
by neighboring electronic states through the Jahn- Teller 
distortion (JTD). It has been concluded that the effect 
of the JTD on the RXS spectra is much larger than that 
of the 3d-4p Coulomb interaction. 

RXS experiments were carried out also on YVO3 1| and 
YTi03^]. Since the JTD is considerably smaller in these 
materials than in LaMnOs, one may think the effect of 
the Coulomb interaction important. This is not the case, 
though; ab initio calculations indicated that the effect 
of lattice distortion on RXS intensities is much larger 
than that of the 3d-4p Coulomb interaction even in these 



materials. 113113 The calculation was carried out within the 
muffin-tin (MT) approximation on the lattice parameters 
determined from the experiment. Since the MT approxi- 
mation averages the potential coming from orbitally po- 
larized 3d states, it works to eliminate the effect of the 
anisotropic terms of the 3d-4p Coulomb interaction on 
the 4p states. We found that the spectra consist of sev- 
eral peaks, to which the GdFe03-type distortion and the 
JTD contribute differently, in good agreement with the 
experiment. 

Recently, an RXS experiment has been carried out on 
KCUF3, in which the intensities on the magiMrtie and 
orbital superlattice spots have been reported.E-3[La The 
magnetic RXS has been studied for typical materials such 
as CoO|n| and NiO[^(J. However, the analyses of_±ha 
experimental data are limited to model calculations.EUEj 
As regards the intensities on the orbital spots, they in- 
crease with going down through the Neel temperature 
in the experiment, and are interpreted as a consequence 
of the increase of-, the 3d orbital order parameter in the 
magnetic phase.EZI This interpretation seems obscure un- 
less it is clarified how the 3d orbital order controls the 
RXS spectra. The purpose of this paper is to eluci- 
date the mechanism of the RXS spectra on the orbital 
and magnetic superlattice spots through ab initio calcu- 
lations. We use the full-potential linearlized augmented 
plane wave (FLAPW) method in the LDA+C7 scheme, 
where the local 3d-3d Coulomb interaction is introduced 
on Cu sites.o 

We obtain the RXS intensity on magnetic spots by 
taking account of the spin-orbit interaction (SOI). The 
calculation has to be accurate, since the intensities are 
about three orders of magnitudes smaller than those on 
orbital spots. This seems the first case to have evalu- 
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ated the magnetic RXS spectra in the ab initio level. We 
reproduce well the experimental spectra as functions of 
photon energy and of azimuthal angle. No intensity ap- 
pears without the SOI, indicating that the intensity arises 
from the orbital polarization, not from the spin polariza- 
tion, in the Ap states. On the other hand, the spectra 
change little with turning off the SOI from the 3d states. 
Since the 3d orbital moment disappears in this condition, 
the 3d orbital moment plays no role to induce the orbital 
polarization in Ap states. This shows a contrast with the 
magnetic RXS in NiO and CoO, where the 3d orbital 
moment plays an important role through the mixing of 
the Ap states with the 3d states of neighboring Cu sites 
and the intra-atomic Ap-3d Coulomb interaction. 

We also calculate the RXS spectra on the orbital su- 
perlattice spats. A similar calculation has already been 
carried out.t§ We come to the same conclusion that the 
RXS intensity is mainly controlled by lattice distortion. 
We add several aspects of making clear the mechanism by 
calculating the spectra with changing magnitude of the 
JTD and in various magnetic phases. We explicitly show 
that the magnetic order itself has little influence on the 
intensity and that the intensity increases with increasing 
JTD strength. In this respect, the increase of the inten- 
sity with passing through the Neel temperature has to 
accompany a considerable increase of the JTD strength 
in the magnetic phase. We also predict a new azimuthal 
angle dependence different from the previously assumed 
one but consistent with the experiment .Ea 

In Sec. II, we study electronic structures in various 
phases with the FLAPW method in the LDA+C/ scheme. 
In Sec III, the absorption spectra and the RXS spectra 
are presented in comparison with the experiment. Sec- 
tion IV is devoted to concluding remarks. 



II. ELECTRONIC STRUCTURES 

KCUF3 has a typical perovskite structure, and under- 
goes a structural phase transition at T w 800 K, be- 
low which the JTD and the associated antiferro-orbital 
(AFO) long-range order emerges with the ordering vec- 
tor Qorb = (tt/g, 7t/6, 7r/c). Figure [l] shows the schematic 
view of the ordering pattern. It is tetragonal with 
a = b = 4. 141 A, c = 3.924A at room temperature; 
three different bond lcngthsi-are given by i = 2.253A, 
s = 1.888A, m = 1.962AMy Since the lattice param- 
eters are not known at low temperatures, we use these 
values in the following calculation. There are two in- 
equivalent Cu sites in a unit cell, which will be called as 
A and B sites. With further decreasing temperatures, 
an antiferromagnetic (AFM) long-range order develops 
below 38 K with the ordering vector Q mag = (0, 0,7r/c). 
Therefore four inequivalent Cu sites, A + , A_, B + , B_, 
exist. The local magnetic moment is pointing to either 
the (110) or (110) direction with the value 0.48 /x B - The 
recent magnetic scattering experiment gives the ratio of 
the orbital angular momentum (L) to the spin angular 
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FIG. 1: Sketch of the ordering pattern of KC11F3. Orbitals 
ip x 2_ z 2 and %f) y 2_ z 2 are antiferro-orbitally ordered with the 
ordering vector Q or b = (n /a,n /b,n /c). Arrows represent 
the direction of magnetic moments with the ordering vector 
Qmag = (0,0, 7r/c). Small solid spheres represent F atoms. K 
atoms are omitted in the figure. 

momentum (S) as L/S= 0.29.0 

We carry out the band structure calculation in the 
LDA+U scheme by setting the parameters of the in- 
traatomic Coulomb interaction between 3d orbitals such 
that U = 3.0 eV and J — 1.0 cV. Wc take account 
of the SOI to deal with the orbital moment, which is 
important for the RXS spectra on the magnetic spots. 
First assuming the AFM long-range order, we obtain an 
orbital long-range order with the ordering vector Qorb! 
the hole numbers inside the Cu MT sphere in the 3d 
states Lp x 2_ z 2 (ip y 2_ z 2) and (p 3y 2_ r 2 (( / 9 3a ,2_ J ,2) are given 
by n x 2_ z 2(n y 2_ z 2) — 0.65 n^ y 2_ r 2(n^ x 2_ r 2) = 0.08 in A± 
(B±) sites, respectively. Note that the difference in the 
values between A + and A_ sites and B + and B„ sites 
is negligibly small. The local magnetic moment is ob- 
tained as 0.79^ with S = 0.34, L = 0.11, L/S = 0.33, 
which is somewhat larger than the experimental one. In 
order to see the effect of the magnetic order on orbital 
polarization, we carry out the same LDA+J7 calculation 
in the presence of the ferromagnetic (FM) long-range or- 
der. We obtain the same type of orbital order, in which 
n x 2_ z 2(n y 2_ z 2) = 0.66 n 3y 2 _ r 2(n Zx 2 _ r 2~) = 0.10 in A (B) 
sites, respectively. Since the hole numbers are similar to 
that of the antiferromagnetic state, we conclude that the 
magnetic order has little influence on the orbital order 
parameter. On the other hand, the orbital polarization is 
rather sensitive to the intraatomic Coulomb interaction. 
The simple LDA (U — 0) gives the orbital polarization 
much smaller than the values mentioned above, that is, 
n x 2_ z 2(n y 2_ z 2) = 0.26 and n 3y 2_ r 2(n 3x 2_ r 2) = 0.11 per 
spin in the A(B) sites in the nonmagnetic (NM) phase. 

Figure || shows the calculated density of states (DOS) 
projected on the d and p symmetries of a Cu site, which 
is defined inside the MT sphere, in the AF phase given 
by the LDA+/7 method. The origin of the energy is the 
bottom of the conduction band. The DOS projected onto 
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FIG. 2: Density of states projected on the s, p and d sym- 
metric states inside the Cu muffin-tin sphere, calculated in 
the AFM phase of the LDA+(7 scheme. The origin of energy 
is the bottom of the conduction band. 
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FIG. 3: Absorption coefficient A(ui) as a function of pho- 
ton energy around the K edge of Cu, in comparison with the 
experimental data of fluorescence (Ref. [if]). The core-hole 
energy is adjusted such that the calculated peak position co- 
incides with the experimental one. 



the d symmetric states is almost concentrated in the re- 
gion of energy less than 2 eV. It has a finite gap ~ 0.82 
eV. The left part of DOS is shown in a 20 times mag- 
nified scale. The most part of probability distributes in 
the interstitial region or on neighboring MT spheres, and 
thereby the magnitude of the DOS projected onto the p 
symmetry becomes much smaller than the 3d DOS. We 
refer to the p symmetric states as Ap states, because the 
wave functions are quite close to the atomic Ap wave func- 
tions within the MT sphere. Of course the "4p" states 
lose the atomic character outside the MT sphere, form- 
ing an energy band with its width being as large as 30 
eV. The Ap DOS shows almost no exchange splitting in 
spite of the magnetic ground state, indicating that the 
3d states have little influence on Ap states. 



III. ABSORPTION COEFFICIENT AND RXS 
SPECTRA 



We define the Is— Ap dipole transition density matrix 
T mm> ( e ) assigned to each Cu site as 



r 2 drr l2 dr'[R*{r)rV^4>nM(r)Y 



T yj> (e) 

'mm'\ 1 / ^ I ' "" ' «" 

x [R* u (r')r"P^UnAr')}S(e-e nM ), (3.1) 

where j stands for sublattices A± or B±. Suffixes to, 
m! represent the x, y, and z axes, which are chosen to 
be parallel to the a, b, and c axes, respectively (see Fig. 
0). The 4> n ,k represents the wave function with the band 
index n, wave- vector k and energy e n ^ which is larger 

than the Fermi energy. Operator Vm projects the wave 
function <pn,)s- on the to (to £ x, y, z) component of the 



p symmetric part on the Cu j site. The Rji s represents 
the Cu Is wave function on the Cu j site. We evaluate 
this density matrix using the FLAPW wave function. 



A. Absorption Coefficient 

We first discuss the absorption coefficient A(u>) around 
the K edge. Neglecting the core-hole potential working 
on the Ap states in the final state of the absorption pro- 
cess, we have the expression for A(u>), 



A{u) K^] I r 2 drR* ls {r)rV 

n,k 



5{lo - e„ :k - £is) 



= t£ (u - e u ) + rW (w - e ls ) + t# (w - e u ), (3.2) 

where £i s is the energy of the Is state. The dependence 
on sublattices can be neglected. A(u) is nearly propor- 
tional to the Ap DOS. Figure|| shows the calculated result 
in comparison with the experiment, where e\ s is adjusted 
such that the position of the main peak coincides with the 
experimental one. It is difficult to determine the absolute 
value of eis within the present status of the ab initio cal- 
culation, since several screening processes are neglected. 
The spectral shape is found in good agreement with the 
experimental curve. If the core-hole potential would be 
taken into account, the intensities in low frequency region 
may be further enhanced. 



B. RXS Spectra 

The conventional RXS geometry is shown in Fig. |J; 
photon with frequency to, momentum q^, polarization fi 
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FIG. 4: RXS geometry. Incident photon with wave vector qi 
and polarization a or n is scattered into the state with wave 
vector q/ and polarization a' or 7r' at Bragg angle 0. The 
sample crystal is rotated by azimuthal angle ip around the 
scattering vector G. 



(= a or 7r) is scattered into state with momentum and 
polarization p! (= a' or 7r'). The scattering vector G is 
defined by q/ — q^. In such a situation, we have the RXS 
intensity I(G, u) 



I(G,u>) cx 



(3.3) 



with 



M w (G,w) = 

1 V (g\x m (j)\A)(A\x' m (j)\9) 
- {E A - E g ) + »T 

pin(out) 



exp(— zG ■ Tj). 

(3.4) 



Here E^ K ""* J is the a; component of polarization vector 
of the incident (scattered) photon, and j runs over Cu 
lattice sites. The \g) represents the ground state with 
energy E g . The |A) represents the intermediate state 
with energy E^; it consists of the excited electron on the 
4p states and a hole on the Is state. T describes the 
broadening due to the Is core-hole lifetime. The RXS 
spectra are rather insensitive with varying T values; we 
set r = 1 eV. The dipole operators x a (j) at site j are 
defined as Xi(j) = x, a^O) = y, Xs(J) = z in the co- 
ordinate frame fixed to the crystal axes with the origin 
located at the center of site j. Hereafter G is expressed 
in unit of (2w/2a, 2ir/2b, 2ir/2c). Just like the absorp- 
tion coefficient, we neglect the core-hole potential in the 
intermediate state. Then the amplitude Af mm /(G,w) is 
expressed in terms of the density matrices: 



M, 



(G,u) = 
— %= y^exp(-zG • r,-) / — 



,(e)de 



- e + e ls 



^■(3.5) 
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FIG. 5: Calculated RXS spectra on a magnetic superlattice 
spot (001) in comparison with the experiment (Ref. |l^). The 
core-hole energy is set the same as in the case of absorption 
coefficient (Fig. |h. 



1. on Magnetic Superlattice Spots 

For the magnetic supe rlat tice spots G = (00m) with rn 
being odd integers, Eq.( |3.5| ) contains the following com- 
bination of density matrices: 




b+) 



-(B- 



>(e) 



(3.6) 



Only the off-|diagpnal elements survives with an antisym- 
metric formEJ'E3 They would vanish without the SOX 
indicating that they arise from the orbital pola rization^ 
in the 4p states. From this expression, Eq.(3.5) becomes 



I(G,u) cx 



(C yz +C ZX ) / de 



77(e) 



LO - £ + Els 



ir 



(3.7) 



E° ut El n . Since the polarization 
C zx is factored out, the photon 



with Cij = E° nt E) 
dependent part C y . 
energy dependence is independent of polarization. 

The actual calculation in the LDA+C7 scheme is rather 
heavy, since the number of atoms in the unit cell becomes 
as large as 20 in the AFM phase. Highly accurate cal- 
culations are required, since the intensity is about three 
orders of magnitude smaller than that for the orbital su- 
perlattice spots, as shown below. It is not easy to sum up 
many k-points in evaluating 77(e), though. Fortunately, 
summing 8x8x8 k-points in the first Brillouin zone 
seems sufficient in comparison with the result for sum- 
ming 6x6x6 k-points. Figure [| shows the calculated 
spectra for G = (001) as a function of photon energy in 
comparison with the experiment .113 The intensity in the 
the a — > a' channel does not appear in agreement with 
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FIG. 6: Azimuthal angle dependence of the peak intensity on 
a magnetic superlattice spot (005) and on an orbital super- 
lattice spot (331), in comparison with the experiment (Ref. 

0)- 



the experiment. Its spectral shape consists of several 
peaks. The highest and next highest peaks correspond 
well to peaks a and b, although their energies are slightly 
shifted. Note that the core-hole energy has already been 
adjusted such that the calculated absorption peak coin- 
cides with the experimental one. We have checked that 
the spectra change little with turning off the SOI from 
the 3d states. This suggests that the 3d orbital moment 
has little influence on inducing the Ap orbital polariza- 
tion. 



gle dependence, 



From Eq. (3.7), we derive explicitly the azimuthal an- 



I(G,u) oc for cr-> a' 



oc 



j cos# cos tp\ 2 for (j — > 7r' 5 (3.8) 



where sin6» = 0.088 for the (001) spot and sin6» = 0.439 
for the (005) spot. Here the azimuthal angle ijj — is 
defined such that the (110) axis is contained in the scat- 
tering plane. Figure [6] shows the peak intensity on (005) 
as a function of in comparison with the experiment. 
The dependence agrees with the experiment. 
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FIG. 7: Calculated RXS spectra on an orbital superlattice 
spot (331) in comparison with the experiment. The thick 
solid and broken lines show the spectra on the AFM and FM 
phases in the LiDA+U scheme, respectively. The thin solid 
line shows the spectra on the NM phase in the LDA scheme. 
The dotted line shows the experimental data (Ref.jj8||). The 
core-hole energy is set the same as in the case of absorption 
coefficient. 



2. on Orbital Superlattice Spots 

For the orbital superlattice s pot s G = (Mm) with h, I, 
and m being odd integers, Eq. ( |3.5| ) contains the following 
combination of density matrices: 



(e)-r^(e) 

7 (e) 
a(e) (5(e) | . (3.9) 
-5(e) 



The off-diagonal elements are about two orders of mag- 
nitude smaller than the diagonal terms. Neglecting such 
terms, we obtain 




I(G,to) oc 
with 



^x 2 —y 1 



de- 



a(e) — /3(e) 



lo - e + eis 



iT 



C x 2 



pOUt 7^1) 



(3.10) 



(3.11) 



Just like Eq. (3.7), the polarization dependent part 



C x 2_ y 2 is factored out, so that the photon energy de- 
pendence becomes independent of polarization. 

Figure shows the calculated spectrum for G = (331), 
as a function of photon energy, in comparison with the 
experiment. The calculation is carried out on both the 
AFM and FM phases in the LDA+J7 scheme, and on the 
NM phase in the simple LDA scheme. All three cases of 
calculation give nearly the same result, indicating that 
the 4p states are little influenced by the details of the 
3d states, in agreement of the previous calculation of the 
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same kind. In particular, the fact that the last case 
gives the similar result in spite of considerably smaller 
orbital polarizations gives a strong support to the mecha- 
nism that the 3c? orbital polarization has little relation to 
the RXS intensity. This is consistent with prftvjnus stud- 
ies on other transition-metal compounds ,E3E3o where 
the RXS intensity are mainly controlled by lattice dis- 
tortion. The effect of the SOI is negligibly small for the 
these cases. The present calculation reproduce well the 
spectra consisting of three peaks, although peaks c and 
e are much smaller than the experimental ones. 

The azimuthal angle dep ende nce for G = (331) is ex- 
plicitly evaluated from Eq.(3.11). It is given by 



7(G,w) oc | cos (3 sin 2V> | 2 for a — ► a' 

oc | sin cos /3 cos 2^ + cos 9 sin (3 sin tp\ 2 



for a 



(3.12) 



with sin6> = 0.364, tan/3 = 3y/2c/a = 4.02. The_ az- 
imuthal angle ip — is again defined such that the (110) 
axis is contained in the scattering plane. The intensity 
has a period of 2tt not n with respect to tp. This rather 
complex behavior seems reasonable, since the scattering 
vector G is not in a direction of high symmetry. This is 
obviously different from the form previously assumed in 
the analysis of the experimental data.lla Figure || shows 
the peak intensity as a function of tp, in comparison with 
the experiment. The present form is consistent with the 
experiment. 

So far, all the spectra are calculated with the lattice pa- 
rameters at room temperature mentioned before. If the 
RXS intensity arises from lattice distortion, it increases 
with increasing JTD strength. To confirm this, we calcu- 
late the intensity with changing the JTD strength. Fig- 
ure H shows the intensity of the peak at hcu = 8995 eV as 
a function of £— s (the values of m and £+s are kept to be 
1.962A and 4. 141 A, respectively). The intensity mono- 
tonically increases with increasing values of £ — s, con- 
sistent with the above observation. It has been observed 
in the experiment that the intensity increases by factor 
2 with temperature going through the magnetic phase 
transition temperature, and this behavior has been in- 
terpreted as a consequence of a strong-coupling between 
orbital and spin degrees of freedom.EZI The result that 
the spectra are mainly controlled by lattice distortion in- 
dicates that the lattice distortion becomes large in the 
magnetic phase. The direct determination of lattice pa- 
rameters have not been carried out yet. 



IV. CONCLUDING REMARKS 

We have calculated the RXS spectra around the K 
edge of Cu in KCUF3 on the basis of the ab initio cal- 
culation, that is, the FLAPW method in the LDA+?7 
scheme. The lattice distortion is introduced as inputs of 
calculation. We have obtained finite intensities on mag- 
netic superlattice spots by taking account of the SOI. The 
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FIG. 8: Peak intensity of RXS spectra on (331) with varying 
JTD strength. The values of m and £ + s are fixed at 1.962A 
and 4. 141 A, respectively. The circle indicated by an arrow is 
the result for the actual value of £ — s at room temperature. 



intensities are about three orders of magnitude smaller 
than those on orbital superlattice spots, consistent with 
the experiment. This is the first result to have evaluated 
the magnetic RXS spectra in the ab initio level. Since the 
intensity disappears without the SOI, the present result 
indicates that the spectra arises from the orbital polar- 
ization, not from the spin polarization, in the 4p states. 
On the other hand, the spectra change little with turn- 
ing off the SOI from the 3d states. This indicates that 
the 3d orbital moment has no role to induce the 4p or- 
bital polarization. This situation is different from our 
previous finding in the study of the magnetic RXS in 
CoO and NiO that the 4p orbital polarization is induced 
by the 3d orbital moment through the mixing of the 4p 
states with the 3d states of neighboring C u . at oms and the 
intra-atomic Ap-3d Coulomb interactionOEj This differ- 
ence may come from the fact that the 3d orbital moment 
is about 1 /3 of that in NiO and 1/10 of that in CoO. Fur- 
thermore, the mixing effect may become smaller than in 
NiO and CoO, since neighboring Cu sites are always in- 
tervened by F atoms. Closely related is the phenomenon 
of the magnetic circular dichroism (MCD) in the, K-edge 
absorption in ferromagnetic metals Fe, Co, Ni£3 In this 
case, the 4p orbital polarization was found to be induced 
by the 3d orbital moment throughJdie Ap mixing with 
the 3d states at neighboring sites PIP 8 ! 

We have also calculated the RXS intensities on or- 
bital superlattice spots. The spectra are independent of 
whether the system is in the AFM phase or in the FM 
phase, and of the orbital polarization in the 3d states. 
We have explicitly shown that the spectra are mainly 
controlled by the lattice distortion through the calcu- 
lation with varying JTD strength. The present result 
is consistent with the previous s tu dies of the RXS on 
transition-metal oxides ,^t3llSllall3 but shows contrast 
with the RXS around the Lm edge in the quadrupole 
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ordering phase of the rare-earth compound CeBg, where 
the sp ectra seems directly controlled by the quadrupole 
order.E3E3ka 
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